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Abstract: The large- a; behavior of the physical evolution kernels appearing in the sec- 
ond order evolution equations of the singlet F2 structure function and of the structure 
j> ■ function in ^-exchange DIS is investigated. The validity of a leading logarithmic thresh- 

old resummation, analogous to the one prevailing for the non-singlet physical kernels, is 
(f-) , established, allowing to recover the predictions of Soar et al. for the double-logarithmic 

contributions (ln*(l — x), i = 4, 5, 6) to the four loop splitting functions Pqg\x) and Pgq\x). 
Threshold resummation at the next-to-leading logarithmic level is found however to break 
down in the three loop kernels, except in the "supersymmetric" case Ca = Cp. Assuming 
a full threshold resummation does hold in this case also beyond three loop gives some in- 
formation on the leading and next-to- leading single-logarithmic contributions (ln 4 (l — x), 
^ i = 2,3) to Pqg (x) and Pg q {x). Similar results are obtained for singlet fragmentation 

functions in e + e~ annihilation up to two loop, where a large-x Gribov-Lipatov relation in 
the physical kernels is pointed out. Assuming this relation also holds at three loop, one 
gets predictions for all large-x logarithmic contributions to the three loop timelike splitting 
function Pg q (x), which are related to similar terms in Pq g (x). 
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1. Introduction 

Although threshold resummation of non-singlet structure and fragmentation functions is 
by now well established [1-4] in QCD, there is still not a comparable understanding of 
the similar problem in the singlet case. Motivated by some recent progress [5-7], this 
issue is addressed anew in the present paper. As in previous studies [8-16] of the non- 
singlet case, the present approach, as well as that of [5], focuses on the properties of 
"physical evolution kernels" [17-19]. At the difference of [5] however, which uses a matrix 
kernel requiring the simultaneous consideration of another process (scalar ^-exchange Deep 
Inelastic Scattering (DIS)) along with photon-exchange DIS, the present work adopts a 
more "intrinsic" point of view which decouples the two processes: namely, it deals with 
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the study of the scalar physical evolution kernels (defined in section 2 and 5) occurring 
in the second order evolution equations [17, 19] associated respectively to the singlet F 2 
and F,)) structure functions. The large-TV behavior of the F2 kernels is investigated in 
section 3, where it is pointed out that at large N there is actually only one independent 
singlet physical kernel. In section 4, an obstruction to a standard form [9] of threshold 
resummation in momentum space is observed in the three loop F2 kernel at the next- 
to-leading logarithmic level, which is found to be removed in the "super symmetric" case 

Ca = Cf- Assuming threshold resummation at Ca = Cf does hold beyond three loop, the 

(3) 

ensuing large- N predictions for the four loop off-diagonal anomalous dimension Jqg (N) 
are derived. The corresponding predictions of [5] for the double logarithmic terms (now 
confirmed in [6,7]) are recovered, and some additional information is obtained concerning 
the single logarithmic terms. A similar study of the F^ kernels is performed in section 
5, yielding analogous large- N information on jgq (N). Fragmentation functions in e + e~ 
annihilation are dealt with in section 6, where a large- x Gribov-Lipatov relation between 
the spacelike and timelike physical kernels is pointed out at two loop. Assuming a similar 
relation is valid at three loop, all large-x logarithmic contributions to the three loop timelike 
splitting function P gq (x) are predicted. The conclusions, and some additional discussion, 
are presented in section 7. The connection between the present approach and that of [5] is 
explained in Appendix A. 

2. Physical evolution kernels of the F 2 singlet structure function 

Consider the singlet F2 structure function in Mellin moment space: 



F(N , Q 2 ) =< e\ > (c q (N, Q 2 /fi 2 , a 9 )q(N, fi 2 ) + C g (N, Q 2 //A a 8 )g(N, /i 2 )) , (2.2) 



where C a (N,Q 2 / /j, 2 ,a s ) (a = q,g) are the singlet quark and gluon coefficient functions, 
q(N,fi 2 ) = YliliiQi + Qi) the singlet quark distribution, g(N,fi 2 ) the gluon distribution, 
and the factorization scale in the coefficient functions has been chosen to be equal to the 
renormalization scale /i 2 with a s = a s (ii 2 )/4iT. Taking two derivatives of eq.(2.2) with 
respect to Q 2 , and eliminating the quark and gluon distributions, one obtains [17,19] the 
second order physical evolution equation: 



J(N,Q 2 ). K and J are renormalization group (and scheme) invariant quantities, which 
are obtained as combinations of coefficient functions. One gets: 




(2.1) 



At large Q 2 , it satisfies the standard OPE representation: 
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R(N q2) = C g (N, Q 2 /» 2 , a s )C g (N, Q 2 /» 2 , a s ) - C g (N, Q 2 / y 2 , a s )C q (N, Q 2 /» 2 , g s ) 

C q (N,Q 2 / f , 2 ,a s )C g (N,Q 2 / f i 2 ,a s )-C g (N,Q 2 / f , 2 ,a s )C q (N,Q 2 / f , 2 ,a s ) ' 

(2.4) 

and 



J(N q2) = CgjN, Q 2 /fi 2 , a s )C g (N, Q 2 /n 2 , a s ) - C g (N, Q 2 /y 2 , a s )C q (N, Q 2 /fi 2 , a s ) 
C q (N,Q 2 /^,a s )C g (N,Q 2 /^,a s )-C g (N 1 Q 2 /n 2 1 a s )C q (N 1 Q 2 / f i 2 ,a s ) ' 

(2-5) 

where C a = dC a /d\n.Q 2 (at fixed /j, 2 ^ Q 2 ). Moreover, the renormalization group equations 
for the coefficient functions yield: 

C q (N, Q 2 /fi 2 , a s ) = P(a s ) dCq{N, ^ /fl2,as) - lqq (N, a s )C q (N, Q 2 /fi 2 , a s )- lgq (N, a s )C g (N, Q 2 /^ 2 , a s ) , 

(2.6) 

dC (N O 2 la 2 a ) 

C g (N, Q 2 //A a s ) = gV ^ ' ^( N , a s )C q (N, Q 2 /fi 2 , a s )- lgg (N, a s )C g (N, Q 2 / f i 2 ,a s ) , 

(2.7) 

where the moment space anomalous dimensions are related to the splitting functions by 
the standard convention: 

Jab(N,a s ) = - f dx x N - 1 P ab (x,a s ) , (2.8) 
Jo 

so that 

sq{N,n 2 )\ = _n qq (N,a s ) lqg {N,a s )\(q(N,n 2 )\ 

\g(N^ 2 y ' \l gq (N,a s ) 7gg (N,a s )J\g(N, f i 2 )J ' 1 ' ) 

Taking one Q 2 derivative of eqs.(2.6) and (2.7) one also gets: 

C q (N, Q 2 /fi 2 , a 8 ) = (3(a s ) dCq{N, ^ /fi2,as) - lqq (N, a s )C q (N, Q 2 /fi 2 , a s )- lgq (N, a s )C g {N, Q 2 /fi 2 , a s ) , 

(2.10) 

and: 



C g (N, Q 2 //A a a ) = P{a s ) dC9{N 'f a ^ 2 ' as) - lqg {N, a s )C q (N, Q 2 /fi 2 , a s )- 7gg (N, a s )C g {N, Q 2 /f?,a s ) . 

(2.11) 

Inserting eqs.((2.6)-(2.7)) and ((2.10)-(2.11)) into eqs.(2.4) and (2.5), then setting fi 2 = Q 2 
(which is legitimate since K and J are renormalization group invariant quantities), and 
using the expansions of the quark and gluon coefficient functions: 
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i=0 
oo 



7ab (iV,a s ) = ^ 7 S(iV)< +1 (2.15) 

i=0 



C q (N,l,a s ) = l + J2 c ?( N )< ( 2 - 12 ) 

i=i 

oo 

C g (N,l,a s ) = J2^\N)ai , (2.13) 

i=i 

as well as those of the beta function and of the anomalous dimensions: 

oo 

/3(a s ) = -J>< +2 (2.14) 

i=0 

with (a, b) = (q,g), one obtains: 

oo 

K(N, Q 2 ) = Y, K{i) ( N H +1 ( 2 - 16 ) 

i=0 

and 

oo 

J(N,Q 2 ) = Y,J (l) (Ny+ 2 , (2.17) 

i=0 

where a s = a s (Q 2 ), and K( l \N) and jW(iV) are expressed as (rather long) combinations 
of the coefficient functions and anomalous dimensions expansion coefficients. 

3. Large-N 

Let us now consider the large- N limit. Using the known large- A expansions of the anoma- 
lous dimensions and of the coefficient functions, one deduces the large-A expansions of 
K®(N) and jW(JV). In the following, I shall focus on K®(N) only, since one can show 
(Appendix A.l) that at large N, J(N, Q 2 ) can be expressed in term of K(N, Q 2 ) and the 
physical evolution kernel K ns (N,Q 2 ) of the non-singlet F2 structure function F ns '. 



J(N, Q 2 ) ~ K ns (A , Q 2 ) - K ns (N, Q 2 ) [K(N, Q 2 ) - K ns (N, Q 2 )\ + 0{l/N 2 ) , (3.1) 
with 

F ns (N, Q 2 ) = K ns (N, Q 2 )F ns (N, Q 2 ) . (3.2) 

Thus at large A, only K®(N) and K$(N) (with K ns {N,Q 2 ) = T,Zo K ™( N ) a i +1 ) are 
independent functions. In particular, eq.(2.3) can be solved at large A as: 



8(F - K ns F) 



, ~K ps {N,Q 2 ){F-K ns F) 



dlnQ 

where K ns F = -jf^F from eq.(3.2), and I defined: 



(3.3) 
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K ps (N, Q 2 ) = K(N, Q 2 ) - K ns (N, Q 2 ) . (3.4) 

K ps can be viewed as a "pure singlet" contribution to K. Since threshold resummation [1,2] 
is well established [9,20] for K ns (N,Q 2 ), one can consider only 1 K(N,Q 2 ). 
a) One loop : 

For the one loop physical kernel coefficient (TV) one finds: 

tf-<°> (AO = - 7 <°> (N) - 7 <«J> (N) - (3 , (3.5) 

which depends on the one loop diagonal anomalous dimensions (and /?o = — | n /)- 

Using the large- A" asymptotics [21]: 

jW(N)~A? +1 \nN-B? +1 (3.6) 
(a = q,g) where In N = In N + 7#, one deduces: 

(AT) ^ K u InN + K w , (3.7) 

with: 

K n = -(Af + Al) (3.8) 

K 10 = B{ + B\ - f3 , (3.9) 

where [22] A\ = 4C F , A\ = 4C A , B\ = 3C F , Bf = fa. One thus gets: 



K u = -4(C A + C F ) (3.10) 
K w = 3C F . (3.11) 

b) Two loop : 

Beyond one loop, the K^(N) 's depend also on the off-diagonal anomalous dimensions, 
as well as on the quark and gluon coefficient functions. Using eq.(3.6) and the large- A 7 " 
asymptotics of the off-diagonal anomalous dimensions [24]: 



^S( N ) ~ ]y£A + i;ln J '# (3.12) 

3=0 
2i 

3=0 

(with D\q = —2nf, Aio = —2C F ), and of the coefficient functions [25]: 

x It should be noted that the large N singlet quark coefficients (eq.(3.14) below) coincide [25] with the 
non-singlet ones, thus threshold resummation of K ns (N,Q 2 ) is related to threshold resummation of the 
singlet quark coefficient function itself. 
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2i 



,(0 



-(0 



i=o 



2i-l 



one gets for i > 1: 



2t 



3=0 



At two loop (i = 1) one finds: 



^22 

K 2 i 



1 i > 

^2i/3o-(^ + ^) 
2cff - ^« (A? - A?) + ±$ (Bf - B\ - /3 0J 



ft 



^20 = tAoA) + + #1) " ft 



2$ + ^$(2*? 



A) , 



with: 



(3.14) 
(3.15) 



(3.16) 



(3.17) 
(3.18) 

(3.19) 



Dji = Dji/rif = -2Dji/D w (3.20) 
$ = = -2cgV^o , (3-21) 

and the brackets contain only one loop quantities. Using known results on one loop coeffi- 
cient functions and two loop anomalous dimensions (see e.g. [23-25] and references therein), 
one gets: 



K 22 = -2{C A + C F )(3 (3.22) 
K 21 = ( - y + 8C2) C\ + ( - y + 8C2) C A C F - y C A f3 - y C F (3 (3.23) 

^20 = (-4 + 12C 3 )Cl + (-10 - 12( 3 )C A C F + (| - 12( 2 + 24Cs)c| 
/27 \ 

+2C A f3 + (y + 10C 2 )C F /? - A , (3.24) 

with /?i = -7C\ - \\C A C F + (5C A + 3C F )/?o. 

c) Three loop : 
At three loop, one finds: 
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K%i = Dsi/3o + (lower loop quantities) (2 < i < 4) (3.25) 
K31 = -D3i/3o — + ^3) + (lower loop quantities) (3.26) 
K 30 = £>3o/3o + (-Bf + Bl) - p 2 + (lower loop quantities) . (3.27) 

Using results on two loop coefficient functions (see e.g. [5, 25, 26] and references therein) 
and three loop anomalous dimensions [23,24], one gets: 



K u = (3.28) 
K 33 = ~(C A + C F )Pl (3.29) 

K 32 = (-^ + 24C 2 ) C 2 A f3 + - 8C2) C A C F (3 - 16C 2 /? - ™ C A (3 2 - ^C F f3 2 

-2(C A + C F )f3 1 (3.30) 

X31 = ( - ^ + ^C 2 - 24C 3 ) C 2 A f3 + ( - ^ + y C2 + I68C3) C A C F /3 + (-64 + 72( 2 - 96C 3 )C 2 & 

/ 238 \ 9 / 313 \ 9 
+ ( - — + 8C2) Ca/?o + ( - -g- + 4C 2 ) Cf0o - 4C A Pi - 5C F /3i - 2/3 /?i (3.31) 

/224 64 > 176 9 \ ~ /2039 64 A 176 ~\ » „ /605 \ _ o 

+ ("g- + yC 2 - — C 2 2 + 176C 3 )ci + (— + yC 2 - —C,1)C\C F +{—- 176( 3 )C A C F . 

d) Four loop : 

For future use, I also quote the large- N structures at the four loop level. One finds: 

3 ~ 

= —D&fio + (lower loop quantities) (2 < i < 6) , (3.32) 

where I restricted to the range 2 < i < 6 which will be of interest in the following. The 
"lower loop quantities" in eqs.(3.25), (3.26), (3.27) and (3.32) are given by very long ex- 
pressions, which need not be written down explicitly here, and depend both on anomalous 
dimensions and (quark and gluon) coefficient functions. One should mention however 
that these expressions do not involve the A^'s coefficients, therefore the K^s (hence 
the K^(Nys at large N) are independent of the off-diagonal j gq (N) anomalous dimen- 
sion 2 . Moreover, the dependence upon the one loop coefficient D\q = — 2nj reveals itself 

(i) 

only through the appearance of the "renormalized" off-diagonal coefficients Dij and c g - 
(eqs.(3.20) and (3.21)). 

4. Momentum space results and threshold resummation 

4.1 An obstruction to threshold resummation 

In this section, a few observations are presented. 

1) The momentum space kernel K(x,Q 2 ) is related in the standard way to the Mellin 
moment space one by: 



2 This is not the case for the K (i) (N)'s at finite N for i > 2. 
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K(N,Q 2 ) = f 1 dxx N ~ 1 K(x,Q 2 ) . (4.1) 
Jo 



Usually, the moment space kernel K(N, Q 2 ) cannot be inverted analytically to momentum 
space (due to the presence of TV-dependent functions in the denominators of eq.(2.4)). 
However, this is no more the case at large N, as is clear from the expressions in the 
previous section. 

2) Using the large-x-large-iV correspondence: 

!^ < _ > L1Z^ 1 „.. +1;V + ..., (4 . 2) 

(l-x)+ n + 1 

the momentum space physical kernel expansion coefficients (x) (obtained by inversion 
of the moment space ones) are seen for < i < 2 to have only a single-logarithmic 
enhancement at large-x (cf. in particular eq.(3.28)), similarly to the non-singlet kernel 
[12,13], namely: 

K „ )(x) _ o( !^£)) (43) 

3) Consequently, it makes sense to ask whether an x space threshold resummation, 
similar to the one valid [9,12] for the non-singlet physical kernel K ns (x, Q 2 ), is also possible 
here, namely whether for x — s> 1, (1 — x)K(x,Q 2 ) effectively depends only on the single 
scale (1 — x)Q 2 : 

where J (Q 2 ) is a renormalization group invariant quantity, so that: 



J ((1 - x)Q 2 ) = ji a s + a 2 s [-j^ L x + j 2 ] (4.5) 
+ allhPlLl - {nfa + 2j 2 p )L x + h] 

+ a A a [-hftLl + (|ji/3i/3 + 3j 2 /3 2 )4 - fafo + 2j 2 /3i + Sj 3 f3 )L x + j 4 ] + ... , 

where a s = a s (Q 2 ) and L x = ln(l — x). Comparing eqs.(4.4) and (4.5) with the exact 
result (known up to 0(a^)), inverting to momentum space the results of section 3 to get 
the leading 0(1/(1 — x)) term in K(x,Q 2 ), one finds : 

i) The leading logarithms in the exact result agree with the leading logarithms in 
eq.(4.5), with 

ji = -K u = A{C A + C F ) . (4.6) 

ii) There is a discrepancy at the next-to- leading logarithmic level, starting at 0(<2g), 
namely: 

Comparing the 0(a 2 ) term in eq.(4.5) with the corresponding term in the exact result, one 
determines 



-8- 



32 = -K 2 i = (y - 8( 2 )(C A + C F )C A + —C A Po + —C F (3 . (4.7) 

However, comparing with the exact coefficient of the single !siL_ term occurring at 0(a? s ) 
in K(x,Q 2 ), one finds the latter is not given by —(jifii + 2j 2 f3o) as required by eq.(4.4) 
and (4.5), but rather by: 

2K 32 = -[j 1 /3 1 + 2(j 2 + 5)p ] , (4.8) 

with 



5 = 16(C A - C F )[(2 - ( 2 )C A - C F ] . (4.9) 

This discrepancy represents an obstruction to a standard threshold resummation at the 
next-to-leading logarithmic level of the singlet physical kernel. The situation here looks 
similar to the one occuring [11, 14] in the case of the (non-singlet) F F structure function 
[27, 28]. 

4.2 A conjectured threshold resummation for C A = C F and four loop predictions 

It is remarkable that the obstruction 8 to threshold resummation at the next-to-leading 
logarithmic level vanishes in the "super symmetric" case C A = C f . This fact suggests that 
a standard threshold resummation, similar to the one available for the non-singlet physical 
kernel K ns (x, Q 2 ), might indeed also be possible for the singlet kernel in this case. 

Let us derive the four loop predictions following by assuming the validity of the ansatz 
(4.4) and (4.5). As a warm-up, I consider first the two and three loop predictions. Eqs.(4.4) 
and (4.5) imply, after converting to Mellin moment space: 

a) Two loop prediction : K 22 = ^ifn/fo, from the 0{a 2 s L x ) leading logarithm (LL) term in 
eq.(4.5), in agreement with the exact result (eq.(3.22)), as expected from the observations in 
section 4.1. Using eq.(3.17), this constraint determines the two loop off-diagonal anomalous 
dimension coefficient D 22 from one loop quantities. One gets: 

D 22 = in f C AF , (4.10) 

where Caf = C A —C F , which agrees as expected (after straightforward inversion to x-space, 
eq.(2.8)) with the exact result [24] for the double logarithmic term in Pqg^ (x). 

b) Three loop predictions : 

i) Km = 0, from the absence of an 0{a? s L^) double logarithmic term in eq.(4.5). Using 
eq.(3.25), this constraint determines D34 from lower loop quantities. One gets: 

D 34 = -^n f C 2 AF , (4.11) 

which indeed agrees as expected with the exact result [24]. 

ii) A 33 = ^KnP 2 , from the 0(a z s L 2 x ) LL term in eq.(4.5). Then eq.(3.25) yields: 

£33 = ln f C AF (3C F - p ) , (4.12) 
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which again agrees as expected with the exact result [24]. Thus double logarithmic terms 

(2) 

(ln*(l — x), i = 3,4) in the off-diagonal three loop splitting function Pg g (x) are correctly 
reproduced. 

iii) K32 = ^(Kiifti + 2K 2 if3 Q ), from the 0(a^L x ) next-to-leading (NLL) logarithm in 
eq.(4.5). Eq.(3.25) then yields: 

D 32 \predict = ^n f C AF [5po + (28 - 24( 2 )C A + 9C F ] . (4.13) 
This time, the above prediction disagrees for Ca / Cp with the exact result [24]: 

D32 = |n/CAF[5A) + (4 - 12Q 2 )C A + 21C F ] , (4.14) 

as expected from the fact that S 7^ 0. I note however that for Ca = Cf, D32 is correctly (as 
again expected) predicted to vanish. Since D34 and D33 also vanish in this limit, it follows 
after inversion to x-space that the leading single logarithmic term (ln 2 (l — x)) in Pgg\x) 
also vanishes for Ca = Cf- This (a priori surprising) fact finds a nice interpretation in the 
present framework: it is seen to follow from the correctness (up to 0(a%)) of eq.(4.4) and 
(4.5) for C A = C F - 

Moreover, the large /?o term in eq.(4.13) is also correct. This result should be expected 
since the discrepancy 5 (eq.(4.9)) is non-leading at large (3q. Thus one might conjecture 
that the ansatz (4.4) and (4.5) is valid at large /3o even for Ca ^ Cf- 
c) Four loop predictions : the ansatz (4.4) and (4.5) similarly implies at the 0(a A s ) level, 
after converting to Mellin moment space: 

i) -K46 = -^45 = 0, from the absence of O(afL^) and 0{a^L^) double logarithmic 
terms in eq.(4.5). Using eq.(3.32), these constraints determines the four loop off-diagonal 
anomalous dimension coefficients D^q and .D45 from lower loop quantities. Using known 
results for the three loop coefficient functions [5,25], one gets: 



D m = (4.15) 
Dtt=nfC 2 AF {-\cF + \fo) . (4.16) 

ii) K 44 = \K lx ffi, from the 0(a^L 3 x ) LL in eq.(4.5). Then eq.(3.32) yields: 

' ' ' (4-17) 
The predictions eqs.(4.15)-(4.17), first obtained in [5], have recently been shown to be 
correct in [6,7]. Thus, after inversion to x-space, one finds that large x double logarithmic 
terms (ln l (l— x), i = 4, 5, 6) in the off-diagonal Pqg\x) splitting function are again correctly 
predicted by the present ansatz. The close connection with the approach of [5] is explained 
in Appendix A. 2. 

iii) K 43 \predict = §Kn/?i/3 + Knffi, from the 0(a 4 s L 2 x ) NLL in eq.(4.5). Eq.(3.32) 
then yields: 
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D 43 \predict = n f C AF [ - ^/3 2 + o Pi(C A , C F ) + P 2 (C A , C F )] . (4.18) 

where Pi and P 2 are respectively linear and quadratic homogeneous polynomials in (C A , C F ), 
which need not be written down since the corresponding terms are not expected to be cor- 
rectly reproduced anyway by the present ansatz. However, the vanishing for C A = C F of 
-D43 is again expected to be a valid prediction. Since D45 and -D44 also vanish in this limit, 
the vanishing of D43 implies, after inversion to x-space, the vanishing for C A = C F of the 
leading single logarithmic term (ln 3 (l — x)) in Pqg\x) (similarly to the three loop case). 
Moreover, the large 0q term should also be correct. 

iv) K 42 \predict = \K XX 2 + K21P1 + 3.K31A). from the 0{a A s L x ) NNLL in eq.(4.5). 
Eq.(3.32) then yields: 

D 42 \predict = n f [Q - \c 2 )c A 2 + ( - p - ^ 2 )c F /3 2 + (26 + 2( 2 - 16C 3 )C F f3 

+ (l8 - 12(2 - ^Cl " 24Cs + 384C 4 )C$ + Caf(/3 Qi(Ca, C f ) + Q 2 (C A , C F ) 
2 81 

+ ~n f f5 , (4.19) 

where Q\ and Q 2 are respectively linear and quadratic homogeneous polynomials in (C A , C F ) 
which again are not expected to be correctly predicted by the present ansatz. The last 
term in eq.(4.19), which depends on 5 (eq.(4.9)), is also obviously incorrect, and should be 
absent in the exact answer, since it is not even a polynomial in the color factors. However, 
this term vanishes for C A = C F , which represents a non-trivial consistency check of the 
present proposal. 

It should be noted that the single logarithm coefficients D 2 ±, D31 and D41 are not pre- 
dicted in the present approach: rather, they represent input parameters for the threshold 
resummation, since they determine respectively K 2 \, K31 and K41, hence the O(o? s ), O{o? s ) 
and 0(cLg) constant terms j 2 (see eq.(4.7)), j$ and j± in eq.(4.5). It is also remarkable 
the leading single logarithmic coefficients D32 and .D43 are both predicted to vanish when 
C A = C F , which may be the signal of a systematic structure. 

5. (ft- exchange DIS 

In the case of ^-exchange DIS, quite similar formulas apply. Defining the moment space 
structure function by: 

F (j> (N,Q 2 )= [ dx x N - 1 F (j> (x,Q 2 ) , (5.1) 
Jo 

we have the large-Q 2 OPE representation [5]: 

F^N, Q 2 ) = C M (N, Q 2 /fi 2 ,a s )A q , nud (N, ^ 2 ) + C^N, Q 2 /fi 2 , a s )A g , nucl (N, fi 2 ) , (5.2) 
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where C^ (iV ,Q 2 / fx 2 ,a s ) (a = q,g) are the (/(-exchange quark and gluon coefficient func- 
tions, A q:fluc i(N, fi 2 ) and A ginuc i(N, /i 2 ) the singiet quark and gluon nucleon matrix ele- 
ments, and a s = a s (/x 2 ). The physical evolutions kernels for this process are defined by: 

F^N, Q 2 ) = K^N, Q 2 )F^N, Q 2 ) + J^(N, Q 2 )F^N, Q 2 ) , (5.3) 
and can be expanded as: 

oo 

K 4> (N,Q 2 ) = Y,Kf(N)a\ +1 (5.4) 

i=0 

and 

oo 

J 4> {N,Q 2 ) = Y,jf{N)a i + 2 , (5.5) 

i=0 

where a s = a s (Q 2 ). As in the photon-exchange case, I shall concentrate in the following on 
K ( p(N,Q 2 ), as J^N^Q 2 ) can be expressed at large N in term of K^{N,Q 2 ) and a large- 
N non-singlet "gluonic" kernel K^^ ns (N,Q 2 ) which satisfies a standard form of threshold 
resummation similarly to the non-singlet quark kernel K ns (N,Q 2 ) (see Appendix A.l). 

5.1 Large-iV 

The main formal difference between the photon and the ^-exchange cases is that the leading 
0(aP s ) (resp. 0(a s )) behavior of the quark (resp. gluon) coefficient functions (eqs.(2.12), 
(2.13)) are interchanged, namely we have: 

oo 

C^(iV,l,a s ) = l + ^c«(Ar)a l s (5.6) 
i=i 

oo 

C M (N,l,a s ) = ^cf q (Ny s . (5.7) 
i=i 

Consequently, at large N the K^(Nys coefficients shall not depend upon the off-diagonal 
7 gff (iV) anomalous dimension, but rather on the 7 sg (iV) one. 
a) One loop : 

At this order, the results are the same as in the photon-exchange case, namely: 

K^W-Ki^N + Kto, (5.8) 

with 



Kf, = K n = -4(C A + C F ) 
Kf = K 10 = 3C F . 



(5.9) 
(5.10) 



b) Two loop : 
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Setting [5]: 



2i 



3=0 



2j-l 



one gets for i > 1: 



2i 



<V)~E<i, lnJAr - 

3=0 



At two loop one finds: 



K. 



1 r 



22 



A 22 /3 



2c 



(i) 



7,2 2 



= -A 21 p -{Al + A\) 



2 ^ 



= -A 20 /3 + (B| + B|) - ft 



2 0?> : 



where: 



.7'' 



A 



-2Aji/A 10 
-2c2,/A 10 • 



(5.11) 
(5.12) 



(5.13) 

(5.14) 
(5.15) 

(5.16) 



(5.17) 
(5.18) 



I note that the above results are obtained from those of section 3 simply by interchanging 
the "renormalized" off-diagonal parameters, i.e. Dji's with the Ajj's, as well as the ("renor- 
malized") photon-exchange gluon coefficients with the ("renormalized") ^-exchange quark 
ones, and also the diagonal photon-exchange quark coefficients with the ^-exchange gluon 
ones, consistently with the remarks at the beginning of section 5.1. Using known results on 
one loop coefficient functions and two loop anomalous dimensions (see e.g. [5, 23, 24] and 
references therein), one obtains: 



•"■22 — 



K. 



-2{C A + C F )fo 
16 



(5.19) 



21 



K. 



20 



( - y + 8C2) C\ + ( - y + 8C 2 ) C A C F - ^C A f3 - y C F /3 - 2/3 2 (5.20) 
(-4 + 12C 3 )Cl + (-10 - 12( 3 )C A C F + (| - 12C 2 + 24C 3 )C| 
+ ( - y + 4C 2 ) C A /3 + (l7 + 6C2) C F p - y pi - fix • (5.21) 
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c) Three loop : 
At three loop, one gets: 



= A3j/?o + (lower loop quantities) (2 < i < 4) (5.22) 
K^-y = A3i/3o — (vlf + A\) + (lower loop quantities) (5.23) 
K$ = A 30 /3 + (-Bf + Bl) - /3 2 + (/ower Zoop quantities) . (5.24) 

Using results on two loop coefficient functions [5] and three loop anomalous dimensions 
[23,24], one finds: 

K| 4 = (5.25) 

Kls = -^(C A + C F )(3 2 (5.26) 

Kf 2 = (y - 24C 2 )ci/3 + ( - ^ + T2^C A C F h + (20 - 32C 2 )C 2 /3 - ^C A f3 2 - yC F /3 2 

-2/3g - 2(C A + Cf)0i (5.27) 
K& = ( - ^ + yC 2 - 96C3)ci/3 + ( " ^ " f C2 + 192C 3 )cUC F/ 3 + (-6 + 48( 2 - ±8( 3 )C 2 F /3 

+ + 36C 2 ) C A /3 2 + ( - ^ - 12C 2 ) C F /3 2 + 6Ca(3\ - 9C F ft - ifofc - y /3 3 (5.28) 

+ (x + > ~ ^ + 176 ^ + (T + yC 2 - l f ti)c\c F + (™ - „,)r/- . 

d) Four loop : 
Finally, I quote the four loop result: 

if^ = -A4j/?o + (lower loop quantities) (2 < i < 6) . (5.29) 

5.2 Threshold resummation and four loop predictions 

Here again, one finds the analogue of the ansatz (4.4), namely: 

f n 2x J± (I 1 ~ X )Q 2 ) ( , on, 

K<t,(x,Q ) ~ '- (5.30) 

(1 - x)+ 

is falsified in an interesting way. Indeed, setting: 

J<t>(Q 2 ) = j<pia s + j^a 2 + j^al + ... , (5.31) 

one determines: 

fa = , (5.32) 

j 02 = —if 21 . (5-33) 



_ 14 _ 



whereas the coefficient of the single ln ^_^ term occurring at 0{a? s ) in K^{x,Q 2 ) is found 
to be given by: 

2^2 = " foift + 2(^ 2 + 5 )/3 o ] , (5.34) 
with a "discrepancy" 5^ with respect to standard threshold resummation: 

J = 4(-5 + 8C 2 )(C A - C F f . (5.35) 

The interesting new feature of eq.(5.35) compared to eq.(4.9) is that the discrepancy is 
now quadratic in Caf = Ca — Cf- As we shall see, this feature implies that the predictions 
following from the ansatz eq.(5.30) are correct at O(o%) not only for Ca = Cf, but also 
up to terms linear in Caf- Indeed one obtains (as in the case of the D^-'s) the correct [24] 
two and three loop predictions: 

A 22 = -ACfCaf (5.36) 

and: 



A 34 = -\c F C\ F (5.37) 
A33 = \c f C A f^ - 15C F + VICaf) - (5.38) 



Moreover one also predicts: 



A 32 \pr edict = C F - 3/3^ + 18C F /3 - 27C F 

62 „ „ / 10 _ . \ „ , 70 



'-foC AF " (y " I6C2) CfCaf -(-§-" 32 C 2 ) Cl F ] , (5.39) 
while the exact result [24] is: 



^3 2 = Cf 



- 3$ + !8C F f3o - 27C 2 F 



-f(3oC AF " (y " I6C2) CfCaf - yCl F ] , (5.40) 

where the mismatch between eq.(5.39) and eq.(5.40) has been highlighted in boldface. 
Indeed it is seen that the mismatch concerns only the term quadratic in Caf- 

Finally, using known results on three loop coefficient functions [5], the four loop pre- 
dictions following from leading logarithmic threshold resummation of K^x^Q 2 ) are: 

A 46 = , (5.41) 
A45 = \c F C 2 AF {^ - 6C F + 8Caf) , (5.42) 
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and: 



A 44 = —CfCaf 



]/3o 2 -y^/3o + fc| + ^/3 C^ + (^-8C 2 )c F C^ + (A + ^ C2 )ci F 

(5-43)' 

which have been first obtained in [5], and shown to be correct in [6,7]. 

Moreover, assuming the validity of the ansatz (5.30) also at the NLL level, one obtains 
the additional prediction: 



C F 

268 
~27 



A 43 \predict = C F - — C F p% + 38C 2 F (3 + ( - + — ^— C2 + — ^— CaJ 



78037 15136 



2048 



27 9 3 

f— - y^)^CfC A f + ( - 8677 + ™C 2 + 2048C 3 )C|C AF 



365 
V~27~ 



+C^ F i?i(/3 , C'a, Cf) 



(5.44) 



where the coefficient R\ of the term quadratic in Caf is a linear homogeneous polynomial 
in (C a ,Cf, A)) which is not expected to be correctly predicted (at the difference of the 
coefficients of the terms linear in CUir). Finally, the ansatz (5.30) at the NNLL level yields a 
prediction for A42. The latter however depends on the three loop non- logarithmic constant 
term A30, which could in principle be derived from the results in [24]. Here I only quote 
one important meaningful feature of the ensuing prediction: 



A 42 \predict = C F P 3 ((3 , C A , C F ) + %C F ^ , 

6 po 



(5.45) 



where P3 is an homogeneous cubic polynomial in (C A , C F , /?o) (as in eq.(5.44)), and the non- 
polynomial contribution is quadratic in Caf (see eq.(5.35)). Eq.(5.45) is thus consistent 
with the assumption that the ansatz (5.30) is correct up to terms linear in Caf- 



6. Singlet fragmentation functions in e + e annihilation 
6.1 Physical kernels 

The (transverse) singlet fragmentation function [3,4,13,19,29-32,36,37,39] F T (N,Q 2 ) 
satisfies a short distance representation analogous to eq.(2.2): 



F T (N, Q 2 ) =< e\ > (C g T (iV, Q 2 /fi 2 , a s )q T (N, ^ 2 ) + Cj(N, Q 2 / f ,a s )g T (N, M 2 )) , (6.1) 

where Cj(N,Q 2 //j 1 2 ,a s ) (a = q,g) are the singlet quark and gluon timelike coefficient 
functions, qT(N,fi 2 ) = Y^i=i(liT + Qir) the singlet quark fragmentation distribution and 
gT(N,fj, 2 ) the gluon fragmentation distribution. The corresponding physical evolution ker- 
nels are defined as in eq.(2.3): 



F T (N, Q 2 ) = K T (N, Q 2 )Ft(N, Q 2 ) + J T (N, Q 2 )F T (N, Q 2 ) 



(6.2) 



In fact, all equations of section 2 take a similar form in the timelike case, except the 
renormalization group equations (2.6)-(2.11), where the indices (qg) and (gq) should be 
permuted, namely: 



q T (N, nJ q (N, a s ) 7 JJN, a A fq T (N, M 2 ) 



g T (N, \rf g (N, a s ) 7 J g (N, a s )J \g T (N, » 



(6.3) 



hence: 



Cj(N,Q 2 /^ 2 ,a s ) = f3(a s ) dC ^ N ^^^ ^ 0>)c ^ Q 2 /^ 2 , a s )-7^(iV, a s )Cj(N, Q 2 /^ 2 , a s ) , 

(6.4) 



Cj(N, Q 2 /fi 2 , a.) = P(a s ) dC ' iN f^ 2,asJ -^ q (N, a s )Cj(N, Q 2 /^ 2 , a s )- 7 J g (N, a s )Cj(N, Q 2 /fi 2 , a.) , 

(6.5) 



and also: 

dC T (N O 2 In 2 a ) 

C%(N, Q 2 /^ 2 , a,) = /3(a a ) qK - * -^ q (N, a s )Cj(N, Q 2 /fi 2 , a s )- 7 l(N, a s )Cj(N, Q 2 /fi 2 , a,) , 

(6.6) 



Cj(N, Q 2 /fi 2 , a s ) = f3(a s ) dC ° iN f* /fl2 ' asJ - 7 J q (N, a s )Cj(N, Q 2 / » 2 ,a s )- 7 J g (N, a s )Cj(N, Q 2 /^ 2 , a 



S) 



(6.7) 



6.2 Large- N Gribov-Lipatov relations 



The timelike coefficient functions are presently known up to two loop [29-32], while the 
timelike singlet anomalous dimensions are known up to two loop for the off-diagonal ele- 
ments [33-35], and up to three loop for the diagonal ones [36,37]. Therefore the physical 
kernels K T and J T can presently be computed only up to two loop [19]. 
The results for the K (N,Q 2 ) kernel at large N can be very simply summarized in the 
form of large-iV Gribov-Lipatov like relations. Setting: 



K T (N,Q 2 ) = J2K^ T (N)a^ , (6.8) 



one finds at large TV: 
a) One loop : 



K (0)T (N) ~ K[ x In N + Kf , (6.9) 

with: 

Kl = K u . (6.10) 
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This result is not surprising, since it is a straightforward consequence (see eqs.(3.8) and 
(3.9)) for a = b (a, b = q,g) of the well-known one loop Gribov-Lipatov relation [38]: 



(0), 



(6.11) 



b) Two loop : More interestingly, although a Gribov-Lipatov like relation is known not 
to be valid at two loop at finite N, neither for the anomalous dimensions [33,39], nor for 
the physical evolution kernels [19], one finds it does hold for the latter at large N. Indeed, 
setting for N — > oo: 



i=0 

one gets, using the results in [29-32] and [33-35]: 

Kl = K 2l (i = l,2) , 



with: 



KL = fa(k> ~ [2c<f - « - A?) 

Kl = lAiiA " (A? + Af 



A) 



2c 



if " l$ T (Af - Af) + Ic<f « - Bf - p ] 



Kl = ± + {Bf + Bf) - 



A) 



~^ T + \c^{Bf-Bf-M 



The A^-'s are defined by: 



with: 



Ag = aJj/Cf = -2Aj-/Af , 



(i > 0), and 



with: 



~(i)T _ (i)T jc _ _ 2 (i)T /a t 

c gj ~ c gj /° F _ 2c si / A io ' 



(6.12) 



(6.13) 



(6.14) 
(6.15) 

(6.16) 



(6.17) 



(6.18) 



(6.19) 



(6.20) 
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(i > 1). Moreover: 



2i 

i=o 

Using eqs.(6.14), (3.17) and (6.15), (3.18), eq.(6.13) is seen to be equivalent to: 



A£ 2 - ^22 = (6.22) 
A21 - ^21 = -(*T5^)(Af - A?) , (6.23) 

where S T sf = f ~ f, and the relations [37,40,41]: S T s^ = 5 T sBf = and [4] &rsc$ = 
SrsCqi = have been used, together with 3 : 

&rs$ = . (6.24) 

One also finds that a relation similar to eq.(6.13) does not hold for the K^q constant term. 
Instead one gets: 

#20 - K 20 = -12( 2 (C A + C F )(3 . (6.25) 

Eq.(6.25) means that in momentum space only the part of K^ T (x) satisfies a 

Gribov-Lipatov relation, not the 5(1 — x) part, namely we have for x — > 1: 

K^(x) ~ ^Hl-x)-K 2X + {K2Q _ C2K22)§{1 _ x) (6 26) 

(1 XJ + 

^ (1)T (x) ~ 2 ^y-^)-^ + (4 - C 2^22)5(l - X) . (6.27) 

Similar large- x relations 4 have been observed [15,16] up to i = 2 for the non singlet physical 
evolution kernels Kn) T and K n % }. 

6.3 Three loop timelike predictions 

Given the non-trivial character of eqs.(6.13) (see in particular eq.(6.23)), it is natural to 
assume they are not accidental, and that similar large- N Gribov-Lipatov like relations also 
hold beyond two loop, as in the non-singlet case [15,16]. At three loop one gets for N — > oo: 

4 

K^ T (N) ~J2 K M ln * N ( 6 - 28 ) 

i=0 

with: 



3 The normalization convention [42] adopted here (see eqs.(6.1), (6.3) and (6.11)) requires the timelike 
gluon coefficient function Cj to be 1/2 the standard one in [29-32] (and the timelike gluon distribution gr 
to be 2n/X the standard one). 

4 Similar relations hold [37,40,41] for the 1/(1 — x) + parts of diagonal splitting functions; in this case, 
the 5(1 — x) terms are also equal. 
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Kji = Ag^o + (lower loop quantities) (2 < i < 4) (6.29) 

Kji = Ag-^o - (Af + Af) + (lower loop quantities) (6.30) 

Kj = A3 /3o + (-Bf T + -£>f T ) — /?2 + (lower loop quantities) , (6.31) 

where the "lower loop quantities" do not depend on the off-diagonal 7^ (A) anomalous 
dimension, i.e. on the Dj^s. Assuming the large- A Gribov-Lipatov relations: 



Kg = K 3t (1 < i < 4) , (6.32) 

eqs.(6.29), (3.25) and (6.30), (3.26) then show that the A^'s (1 < i < 4) can be predicted 
in term of the D^s (and three loop cusp anomalous dimensions) and lower loop quantities. 
Indeed eq.(6.32) yields: 



A34 



D u = 
A33-A53 = (Af 

A32 ~~ ^32 = 



- A\)[- (6 TS c$) + (6 T scW)c$ 



^Al-A^(5 TS ^)(^ + ^) 
+\(A{ - A\)(5tsc^)(K t 21 + D 21 ) + A(5 7 
-(A{ - A\)(5 TS ~ef) + (B{ - B{ - 2A,)(* rfl <§ ) ) 



5tsc ( $) 



-(Alo-^o) 



D 22 + (Af - A\)~c^ 



+ 



(*r*$) 



K - A?)c« - 4cg 
-\(A\ - A\)(5 TS c§)(B{ - B\ - 
-\(Al-Al)(^ T -D 21 c^ 
+(Stsc { £) [Id 22 (B1 - B\ - p ) + (A{ - A\)(B{ - B\) 
+ (A{ - Aftc™ - (Bf - B\ - 2A,)cg - MAf - A\)~c 



,(i) 

'<?1 



;(1) 



(6.33) 
(6.34) 



(6.35) 



(with a similar, even more lengthy equation for A 3^ — D31, not written down for brevity), 
where I used [4] ^tsc^ = &TSC q 2 } = 0, and also eqs.(6.22) and (6.23), together with: 



(6.36) 



Thus the large-A Gribov-Lipatov relations eq.(6.32) determines the A^'s in term of the 
Z?3i's and lower loop quantities. One finds: 
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A^ 4 = --C F C l AF (6.37) 



A33 = 3C 'fC 'af(12C a - 9C F - /3 ) (6.38) 



Ais = \c F C AF [{& - 24( 2 )C A + f>lC F + 13/3 )] (6.39) 

= - 10C 2 )^/3 + ( - ^ + 14< 2 ) C^CWfc + + 46C 2 - 5 -f(s)c A C F 

+ (- Q f- 32C 2 + ^C 3 ) C\C F + ( - 14 - 26 C2 + ^C 3 ) C$ . (6.40) 

Inverting to momentum space, one gets for x — > 1: 

4 

^ )T (aO~£Ai;in*(l-z), (6.41) 



i=0 



with: 



A34 = (6.42) 
A33 = \c F C AF {V2C A - 9C F - /3 ) (6.43) 
Al, = -\c F C AF [(S - 12Q 2 )C A + (51 - 12C 2 )C F + 13/3 )] (6.44) 
&31 = " 12C 2 )C|A) + ( - ^ + 16C 2 )c A C F /3 + + 88C 2 - 208Cs) C A C£ 

+ ( - — - 56C 2 + 8OC3) C\C F + ( - 14 - 44C 2 + 128C 3 ) C F . (6.45) 

Of course, the double logarithmic terms (i = 3, 4) could also have been obtained from 
an alternative ansatz similar to the one used in section 4.2. The present Gribov-Lipatov 
ansatz however yields in addition the single logarithmic terms A^ (i = 1,2). I note the 
predicted A^'s vanish when Ca = C F for i = 2, 3,4, similarly to the -D3i's. 
The Gribov-Lipatov relations eq.(6.32) are not expected to be valid for the constant terms 
(i = 0). Nevertheless, it is interesting to speculate about a simple structure of the difference 
K^ — K30, similar to the one observed in the two loop case, eq.(6.25). Indeed one finds: 



- K 30 = -12( 2 (C A + C F )f3 1 + [Aio - D 30 ]f3 



(6.46) 



+ 



(-6 + 48C 2 + 144Cf - 24C 3 )Ci + (8 + 44( 2 - 240C 2 2 - 24( 3 )C A C F 



+(18 - 186C 2 + 192C! + 48C 3 )C|]/3 + (4 - 12C 2 )C A + (6 - 134C 2 )Cf Pi , 



where all explicitly calculated contributions arise from less then three loop, and D 3 q could 
be in principle extracted from the results in [24]. The first term on the right-hand side 
of eq.(6.46) reveals that indeed the difference Kj[ — K 3 q cannot vanish, since this term 
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cannot be canceled by the other contributions, which are all proportional to However, 
as suggested by the right-hand side of eq.(6.25), a number of cancellations may take place. 
In particular, one might expect that all terms with no factor of £2 cancel on the right-hand 
side of eq.(6.46): 

Kj - K 30 = 0(( 2 ) , (6.47) 

a possibility which would be interesting to check against the forthcoming exact result [42] . 

It is clear that similar methods may be used to predict the large- x logarithmic coeffi- 
cients (1 < i < 4) in the three loop off-diagonal timelike splitting function Pgp T ' (x): 

P£ )T (x)~J2DTW(l-x), (6.48) 

i=0 

assuming the analogue of the large- N Gribov-Lipatov relations eq.(6.32) for the ^-related 
physical kernels: 

Kff = K* (1 < i < 4) . (6.49) 

Eq.(6.49) yields relations between the -D^'s and the A3j's, which however depend on the 
one and two loop 0-decay timelike quark and gluon coefficient functions. The latter could in 
principle be determined [37] by analytic continuation of corresponding spacelike quantities. 

7. Conclusion 

The large- x / large- N behavior of the physical evolution kernels associated to the second or- 
der evolution equations for the singlet F2 (photon-exchange DIS) and (</>-exchange DIS) 
structure functions has been investigated. It was shown that for each process at large N 
there is actually only one independent "scalar" physical kernel K(x,Q 2 ) (resp. K ( p(x,Q 2 )) 
for singlet evolution, in addition to the non-singlet one K ns (x,Q 2 ) (resp. K ( f >!ns (x,Q 2 )). 
The singlet kernel K(x,Q 2 ) was found to satisfy up to three loop a leading logarithmic 
standard form of threshold resummation, analogous to the one valid for the non-singlet 
kernel, and which holds presumably also beyond three loop. This assumption allows to 
predict the double logarithmic contributions to the four loop off-diagonal splitting function 
Pqg\x), thus recovering some of the results of [5]. It was shown that this agreement is not 
accidental, and that the assumption of leading logarithmic threshold resummation of the 
singlet "scalar" kernel K(x,Q 2 ) is closely related to the assumption (now confirmed [6,7]) 
of single logarithmic behavior of the "matrix" physical kernel studied in [5]. The present 
more intrinsic "scalar" approach allows however to go a little further, and provides some 
information on the two leading large x single logarithmic contributions (ln*(l — x), i = 2, 3) 
to Pqg\x). It is found that there is an obstruction to threshold resummation (as opposed to 
the non-singlet case) at the next-to-leading logarithmic level, which first manifests itself in 
the three loop kernel (a similar fact has been observed [11,14] in the case of the non-singlet 
longitudinal Fl structure function [27, 28]). This obstruction turns out to be removed 
in the "supersymmetric" case Ca = Cf, as well as at large /3q, which suggests a full 
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threshold resummation might still be possible in these cases, leading to the above mentioned 

additional predictions. On the other hand, no prediction is available for the C(ln(l — x)) 

(3) 

term in Pq g '(x), which should be viewed as an input for the threshold resummation. Of 
course, more work is needed to justify the assumed threshold resummation for Ca = Cf (for 
which a non-trivial consistency check has been provided, see the comment after eq.(4.19)), 
and to better understand the general Ca ^ Cf situation. 

The analogous study of (^-exchange DIS, where the scalar <ft is directly coupled to gluons, 
gives information on the other off-diagonal splitting function P gq (x). The interesting new 
feature in this case is that the discrepancy at three loop with threshold resummation is 
found to be quadratic in Ca — Cf, allowing to predict more terms in the single logarithmic 

(3) 

contributions to Pg q (x). 

A similar approach has been applied to the study of the large- x behavior of e + e~ fragmen- 
tation functions. A large- a; Gribov-Lipatov relation (similar to the one holding [15, 16] for 
the non-singlet kernels) has been observed in the two loop physical kernels. Assuming a 
similar relation is valid at three loop, all large- x logarithmic contributions to the three loop 
timelike off-diagonal splitting function Pg Q (x) are predicted and related to corresponding 
terms in Pqg (x) , which could be checked against exact calculations [42] in the near future. 

Acknowledgements I thank A. Vogt for a very helpful correspondence on the normal- 
ization conventions for the timelike gluon coefficient function, and A. Mitov for raising a 
thought-stimulating question. 



A. Connection with the matrix physical evolution kernel for the (F 2 ,F ( p) 
system 

An standard alternative method [17] to construct physical evolution kernels for singlet 
evolution consists in the simultaneous introduction of another process in addition to the 
photon exchange one, in particular the exchange of a scalar (ft directly coupled to gluons. 
This system has been considered in [5] . The matrix physical evolution kernel for the coupled 
(F 2 ,F ( p) system is defined by: 

/F(A,Q 2 )x (K 22 (N,Q 2 )K 2(f> (N,Q 2 )\(F(N,Q 2 )^ 

VF (A,Q 2 )/ \K^{N,Q 2 )K H {N,Q 2 y\F^N,Q 2 )) ' K ' } 

with F as in eq.(2.1). Taking the Q 2 derivative of both sides of eq.(A.l), considering the 
first line (relevant for the photon exchange case), and eliminating F^ and Fq, one recovers 
the "scalar" evolution equation eq.(2.3) with the identifications: 

K{N, Q 2 ) = K 22 (N, Q 2 ) + K H {N, Q 2 ) + (A.2) 
J(N,Q 2 ) = K 2( p(N, Q 2 )K^, 2 (N, Q 2 ) - K 22 (N, Q 2 )K^(N, Q 2 ) 
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Taking next the derivative of the second line of eq.(A.l), and eliminating F and F, one 
recovers similarly eq.(5.3), with K<j,(N, Q 2 ) and J<p(N, Q 2 ) obtained by relations analogous 
to eqs.(A.2) and (A. 3) with the indices 2 and <p on the right-hand sides interchanged. 

A.l Relation between the J (resp. J^) and K (resp. K^) singlet physical kernels 

Consider now the large N limit. I first observe that in this limit, J(N, Q 2 ) can be expressed 
in term of K(N, Q 2 ) and K 2 2(N,Q 2 ). Indeed, it is known [5] that at large N, K 22 (N,Q 2 ) 
and K H {N,Q 2 ) are O(N ) (up to logarithms of N), while K 2(j> {N,Q 2 ) and K^N^Q 2 ) 
are 0(1/N). Thus for N — > oo one can neglect the first term on the right-hand side of 
eq.(A.3), and get: 



J(N, Q 2 ) ~ K 22 (N, Q 2 ) - K 22 (N, Q 2 ) K^N, Q 2 ) + 



K 24> (N,Q 2 
K 24> {N,Q* 



Using eq.(A.2), one then finds: 



J(N,Q 2 )~ K 22 (N,Q' 



K 22 (N,Q 2 ) K(N, Q 2 ) 



K 22 (N,Q 2 



(A.4) 



(A.5) 



Moreover, it follows from [5] that at large N, the O(N ) part (modulo logarithms of N) 
of K 22 (N,Q 2 ) coincides with the non singlet kernel K ns (N,Q 2 ), which is known [9,20] to 
satisfy a standard form of threshold resummation. Indeed, in momentum space one can 
check that eqs.(4.13)-(4.16) of [5] can be cast in the form of eqs.(4.4) and (4.5) (truncated 
to ©(of)), namely we have at large x: 



K 22 (x,Q 2 ) 



: Jns((l-x)Q 2 ) . , kJ 
(1 - x)+ 



K ns (x,Q 2 ) 



(A.6) 



where J ns ((1 — x)Q 2 ) satisfies the analogue of eq.(4.5) (with j ns ^ = 4Cf)- Thus we 
obtain: 



J(N, Q 2 ) - K ns (N, Q 2 ) - K ns (N, Q 2 ) K(N, Q 2 ) - K ns (N, Q 2 ) 



(A.7) 



which, for N — > oo, expresses 5 J(N,Q 2 ) in term of K(N,Q 2 ) and K ns (N,Q 2 ) up to 
0(1/N 2 ) corrections. It thus suffices to study the large- N properties of K(N,Q 2 ). 

Similarly, J<f,(N, Q 2 ) can be expressed at large N in term of K<p(N, Q 2 ) and K^{N, Q 2 ): 



J^N, Q 2 ) - K H {N, Q 2 ) - K H {N, Q 2 ) K^N, Q 2 ) - K H {N, Q 2 ) 



(A.8) 



The results of [5] also imply that K^N ,Q 2 ) satisfies at large N a standard non-singlet 
form of threshold resummation 6 , and can be identified to a non-singlet "gluonic" physical 
kernel. Indeed, one can check that eqs. (4.18)-(4.21) of [5] can also be cast in the form of 
eqs.(4.4) and (4.5) (truncated to 0(af)), namely we have: 



5 For K(N,Q 2 ) = K ns (N, Q 2 ), eq.(A.7) becomes J(N,Q 2 ) ~ K ns (N, Q 2 ), which also follows by taking 
the Q 2 derivative of eq.(3.2), and in this case is valid at finite iV too. 

6 This fact is closely related to the known [5,43] soft gluon exponentiation of C^>, 9 . 
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K H {x, Q 2 ) ~ J ^ ns ( [ {l , X)Q2) ~ K^x, Q 2 ) , (A.9) 
(1 — x) + 

where J^ ns ((1 — x)Q 2 ) satisfies the analogue of eq.(4.5) (with jf ,ns = AC a)- 
To summarize, at large N there are only two independent singlet "scalar" physical kernels 
(K(N,Q 2 ) and K^N^Q 2 )), and two independent non-singlet kernels (K ns (N,Q 2 ) and 
Kcp,ns(N,Q 2 )), which match the four elements of the matrix kernel (A.l). 

A. 2 Single logarithmic enhancement of K 2( j > (resp. K^ 2 ) and leading logarithmic 
threshold resummation of K (resp. K^) 

Setting: 



K 24> {N,Q 2 ) = Y,K§(N)a^ , (A.10) 
=o 



let us first assume the large- N behavior of the off-diagonal kernel coefficients K^^N) to 
be doubly logarithmic (as a priori expected from the behavior of the off-diagonal splitting 
functions): 



3=0 

On the other hand, for the diagonal kernels coefficients Kal{N) (a = 2,<p), assuming 
a single logarithmic enhancement, as follows from their large- iV identification with non- 
singlet physical kernels, and their ensuing threshold resummation properties (Appendix 
A.l), one has: 

oo 

K aa (N,Q 2 ) = Y,K^ a {N)a\ +1 , (A.12) 

with: 



i=0 




i+1 






\n j N 


3=0 




i+1 






ln j N 


3=0 





~ > >7VV" J V (A.l:!) 

K®{N)~^K*$h,>N. (A.14) 

3=0 

Using eqs.(A.10)-(A.14) into eq.(A.2), one deduces the large-iV behavior (eqs.(3.7) and 
(3.16)) of the KW(N)'a (eq.(2.16)), namely: 
i) At one loop: 

K n = K™ + Kf{ ns (A.15) 

K 10 = + Kf 6 ns - p . (A.16) 
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ii) At two loop: 



where: 



iii) At three loop: 



K22 


— -"-22 


+ #22™ " 




K21 


— -"-21 


-"-21 


^21 A) 


K 20 


— -"-20 


-"-20 


^OA) 



4 = 4m = 2 4/^ 



10 • 



(A.17) 
(A.18) 
(A.19) 

(A.20) 



Vk + ( 2^22 



1 



33 



+ 



K33 
K 3l 

K 30 = K% + K& na + 
iv) At four loop: 



A) 



1 ^ 



-^33 + 2 ^22-^21 



-"-32 "+" -"32 



+ Kg* + 



^2 + (\Di) 2 + ^ 2 ^o]Ao - ^2 A 



-^30 + 



1 ^ 
2 

1 * 



n'P + _ n<P 
u z\ + 0^21-^20 



1.-.. 
2 J 
1 

2' 



A) ~ ^D^t 



(A.21) 
(A.22) 
(A.23) 
(A.24) 
(A.25) 



K m = 
K 45 = 



2 46 ^ 4-^34-^22 



In 



22 



A 







- Ibt + + -MM, - l&L 



2 45 



^-^34 -^21 ^-^33-^22 2' 



(^ 







if 44 = ^ + ^44 

3 



+ 



2 22 



2 44 ^ ^"^34-^20 ^ ^"^33-^21 ^ ^-^32-^22 

Di+(iD| 2 ) 2 ]ft 



+ 



(A.26) 
(A.27) 

(A.28) 



with similar, but longer relations (not written down for brevity) for {i < 3). 

If one now takes into account the known [5-7] single logarithmic behavior of K^}{N), i.e 



the fact that Df- = for j > i, one gets: 



K 34 = K 45 = K m = 0, 



(A.29) 



and: 



vi-1 



(A.30) 

Using the relations [5] K% s = -AC F (3 l Q - l /i and K% ns = -4C j4 ^" 1 /«, which follow (at 
the leading logarithmic level) from the threshold resummation of the diagonal kernels 
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mentioned in Appendix A.l, one then recovers the results following from the assumed 
leading logarithmic threshold resummation of the singlet kernel K (section 4.2). These 
results show that the leading logarithmic threshold resummation of K follows from the 
single logarithmic enhancement of the off-diagonal physical kernel K 2l j > , together with the 
(leading logarithmic) threshold resummation of the diagonal kernels K22 and K^. I note 
that eqs.(A.15)-(A.25) allow also to compute the K^s, and thus check the results of section 
3, using the results for the Dfj's in [5]. 

Quite similarly, one can show that the leading logarithmic threshold resummation of 
follows from the single logarithmic enhancement of the off-diagonal physical kernel 
(together with the leading logarithmic threshold resummation of the diagonal kernels K22 
and K^). The (now confirmed [6,7]) ensuing predictions (eqs.(4.15)-(4.17)) of the double 
logarithmic off-diagonal four loop coefficients (4 < i < 6), as well as the similar 
predictions (eqs.(5.41)-(5.43)) of the off-diagonal four loop coefficients (4 < i < 6) are 
thus equivalent to those of [5] . 
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